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ABSTRACT 

We investigate the X-ray number counts in the 1-2 Ms Chandra Deep Fields (CDFs) to determine the contri- 
butions of faint X-ray source populations to the extragalactic X-ray background (XRB). X-ray sources were 
separated into Active Galactic Nuclei (AGN), star-forming galaxies, and Galactic stars based primarily on 
X-ray-to-optical flux ratios, optical spectral classifications. X-ray spectra, and intrinsic X-ray luminosities. 
Number-count slopes and normalizations below 2 x 10"'^ erg cm~^ s~' were calculated in each band for 
all source types assuming a single power-law model. We find that AGN continue to dominate the number 
counts in the 0.5-2.0 keV and 2-8 keV bands. At flux limits of « 2.5 x 10"'^ erg cm"^ s"' (0.5-2.0 keV) 
and w 1.4 X 10"'^ erg cm"^ s"' (2-8 keV), the overall AGN source densities are 1166^1^ sources deg"^ and 
4558;|;207 sources deg~^, respectively; these are factors of ^ 10-20 higher than found in the deepest optical spec- 
troscopic surveys. While still a minority, the number counts of star-forming galaxies climb steeply such that 
they eventually achieve source densities of llll^^J^ sources deg"^ (0.5-2.0 keV) and sources deg"^ (2- 

8 keV) at the CDF flux limits. The number of star-forming galaxies will likely overtake the number of AGN at 
' (0.5-2.0 keV) and dominate the overall number counts thereafter. Adopting XRB flux 

ideg-2 



~ 1 X 10 erg cm ^ s 
densities of (7.52±0.35) x lO'^^erg cm"^ s"' deg"^ 



(0.5-2.0 keV) and (2.24 ±0.11) x lO"" erg cm'^ s" 



(2-8 keV), the CDFs resolve a total of 89.5!§ ;'^% and 86.9!^:5% of the extragalactic 0.5-2.0 keV and 2-8 keV 
XRBs, respectively. AGN as a whole contribute «83% and «95% to the these resolved XRB fractions, respec- 
tively, while star-forming galaxies comprise only «3% and «2%, respectively, and Galactic stars comprise the 
remainder. Extrapolation of the number-count slopes can easily account for the entire 0.5-2.0 keV and 2-8 keV 
XRBs to within statistical errors. We additionally examine the X-ray number counts as functions of intrinsic X- 
ray luminosity and absorption, finding that sources with Lo.5-8 keV > lO"'-''^ erg s"' and A^h < 10^^ cm"^ are the 
dominant contributors to the 0.5-2.0 keV XRB flux density, while sources with Lo.5-8 keV = 10"*^^ -10'*'*^ erg s"' 
and a broad range of absorption column densities primarily contribute to the 2-8 keV XRB flux density. This 
trend suggests that even less intrinsically luminous, more highly obscured AGN may dominate the number 
counts at higher energies where the XRB intensity peaks. Finally, we revisit the reported differences between 
the CDF-North and CDF-South number counts, finding that the two fields are consistent with each other except 



for 2-8 keV detected sources below F2-8 keV 
«3.9f7. 

Subject headings: cosmology: observations - 
1. INTRODUCTION 



1 X 10 erg cm s , where deviations gradually increase to 
- galaxies: active — galaxies: starburst — X-rays: galaxies 



Deep Chandra and XMM-Newton observations have now 

resolved the vast maior ity of the X-ray ba ckground (XRB) 

below ft!8 keV (e .g., ICowie et al.l I200I hereafter C02 ; 

iMoretti et all I2003L hereafter M03: ' Worslev et alJ l2()04l) . 

with much of the remaining uncertainty in the resolved 

fraction attributed to deviations in the absolut e value of 

the XRB i tself due to l arge-scale structure (e.g., iGilli et aP 

12003; Yan g et al.l 120031) and significant instrumental cross- 
^^^^^^^ ^^—^^^ — 

calibration uncertainties (e.g.. lDe Luca & Mo lendi 2004. and 
references therein). The number counts in both the 0.5- 
2.0 keV (soft) and 2-8 keV (hard) bands can be fitted 
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with characteristic double power-law shapes with breaks 
around lO-'^-lO"'"* erg cm'^ s"' (e.g., C02; M03). In- 
tensive optical follow-up campaigns have shown that ex- 
tragalactic sources at bright X-ray fluxes (> 10~'-'-10~'^ 
erg cm"^ s~') are generally unobscure d or m ildly ob- 
scured active galactic nuclei (AGN; e.g. , iBade e t al. 19981 
[Schmidt et al." 1998"; 'Akiyamaet al.'"2003"), while below this 
level several other populations emerge such as obscured AGN 
(e.g., lAlexander et al. 2001; Tozzi et al. 2001; B arger et^TI 
120021 Mainieri et al. .2002i) and s tarburst and quiescent 
galaxie s (e.g., Giacc om et alJl200U lAlexander et al.ll2002at 
iBauer et al.i2002ai:iHornschemeier et alJ2003 !). 

This paper builds upon the overall number-count re- 
sults of C02 and M03, which were both based on the 
1 Ms Chandra Deep Field (CDF) datasets, and extends 
the quiescent galax y nuni ber-count results presented by 
iHornschemeier et'al] i200 3') by investigating how various 
source populations contribute to the XRB using the multi- 
wavelength datasets of the 2 Ms CDF-North (CDF-N) and 
1 Ms CDF-South (CDF-S). In addition to deep X-ray ob- 
servations, these legacy fields have deep HST, radio, and 
ground-based optical imaging, as well as several thousand 
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Fig. 1 . — The CDF-N (black curves) and CDF-S (grey curves) solid angle 
coverage for a 2u flux limit. Solid and dashed curves represent solid angle 
coverage in the 0.5-2.0 keV and 2-8 keV bands, respectively. 

spectroscopic redshifts, allowing classification of different 
source types. We describe our X-ray sample in f|2l while 
our method for estimating incompleteness and bias is out- 
lined in Sj3] Finally, we present the selection and classifica- 
tion of sources and number-count results in S|3 Throughout 
this pape n we adopt Hn = 70 km s"' Mpc"\ 57m = 0.3, and 
J7a = 0.7 JSpergel et alJ2003h . Unless explicitly stated other- 
wise, quoted errors are for a Icr (68%) confidence level. 

2. X-RAY SAMPLE 

Our sample is derived from the catalogs of lAlexander et al.l 
J2003L hereafter A03), which consist of 503 X -ray sources m 
the 2 Ms CDF-N and 326 X-ray sources in the 1 Ms CDF-S. 
We have chosen to combine these two samples to achieve a 
more accurate census of the X-ray source population. We rec- 
ognize that the number counts of these two fields have been 
reported to differ by ~ 30% at faint X-ray fluxes (e.g., C02; 
M03), and we briefly investigate such differences in S^l The 
on-axis sensitivity limits for the CDF-N and CDF-S are w 
2.5 X 10"'^ erg cm"^ s"' and w 5.2 x 10"'^ erg cm"^ s"' (soft) 
and « 1 .4 X 10""' erg cm"^ s"' and w 2.8 x lO""' erg cm"^ s"' 
(hard), respectively. X-ray source fluxes are taken directly 
from these catalogs and include corrections for vignetting, in- 
dividual spectral slopes when known (otherwise F = 1 .4 is as- 
sumed), and contamination of the ACIS blocking filters. 

We have additionally corrected the fluxes for the interven- 
ing Galactic column densities.^ For our number-count esti- 
mates in each band, we have imposed a flux cutoff to the 
above sample based on the 2a limiting flux maps for the 
CDF-N and CDF-S, which were constructed following the 
prescription in §4.2 of A03. This significance level empiri- 
cally matches our limiting sensitivity across the field and is 
adopted to remove a small number of sources with large and 
highly uncertain completeness and flux bias corrections (i.e., 
< 10% completeness; see [J3}- Figure ^ presents the sensi- 
tivity and sky coverage of the CDFs. Twenty soft-band and 
10 hard-band sources were rejected because they were be- 
low this flux threshold. In total, we used 724 sources de- 
tected in the soft band and 520 sources detected in the hard 
band to estimate the X-ray number counts, with sky coverages 
ranging from 0.232 deg^ to 0.004 deg^ depending on X-ray 

' The Galactic column densities toward die CDF-N and C DF-S are (1.3 dz 
0.4) X 10^° cm~^ and (8.8 ± 4.0) X lO" cm"^^ respectively ILockmanl200l 
IStarket alim% . 



flux. Note that the CDF sample presented here consists only 
of point sources and does not include a ny obvious contribu - 
tion from X-ray clusters or groups (e.g., iBauer etal ■I2002bh . 
Thus our total X-ray number counts and resolved XRB frac- 
tion estimates will likely underestimate the real quantities by a 
few percent (for e stimates of this small contribution see, e.g., 
iRosati et all2002l) . 

3. CHANDRA SIMULATIONS 

To understand the effects of incompleteness and bias, we 
created 200 Monte Carlo simulated observations in both the 
soft and hard energy bands for each of the CDFs. We began 
by creating template background images for each field in the 
soft and hard bands following §4.2 of A03. To this end, we 
masked out all 829 known point-sources using circular aper- 
tures with radii twice those of the « 90% point spread func- 
tion (PSF) encircled-energy radii. We filled in the masked 
regions for each source with a local background estimated 
by making a Poisson probability distribution of counts using 
an annulus with inner and outer radii of 2 and 4 times the 
«90% PSF encircled-energy radius, respectively (for further 
details see A03). The resultant images include minimal con- 
tributions from detected point sources while still providing re- 
aUstic contributions from extended sources (e.g., Bauer et aQ 
'2002b), which will cause a slight overestimation of the mea- 
sured background close to extended sources. 

To these template background images we added simulated 
sources at random positions. The fluxes of these simulated 
sources were drawn randomly from the total number-count 
models of M03 between lO-^^-lO"" erg cm s (soft) and 
10~'^-10~" erg cm"^ s"' (hard), respectively. These fluxes 
were converted to X-ray count rates assuming a F = 1 .4 power 
law and the maximum effective area of the CDFs. Assuming 
a somewhat different average X-ray spectral slope (F = 1-2) 
does not dramatically change our results. We have not at- 
tempted to account for any scatter in the distribution of F and 
note that individual sources with significantly steeper or flatter 
spectral slopes are likely to have much higher flux thresholds 
for detection due to the energy dependence of Chandra's ef- 
fective area. The number of CDF sources with deviant X-ray 
spectral slopes is small, and therefore are unlikely to have a 
large effect. 

Exposure times for the simulated sources were derived from 
their positions on the CDF exposure maps (see §3.1 of A03) 
and used to convert count rates to counts. To include the ef- 
fects of Eddington bias (i.e., the measured flux is higher than 
the actual flux due to statistical fluctuations), the counts for 
a simulated source were then redrawn from a statistical error 
distribution. Note that we kept track of the counts both be- 
fore and after including Eddington bias in order to decouple 
the effects of our photometry from those of Eddington bias 
(see Fig.|3jl. Finally, counts for each source were drawn ran- 
domly from a PSF probability distribution function to simu- 
late a real source and then added to the template image. To 
mimic the complex PSF of the multi-observation CDFs, we 
adopted the combined model PSF from the nearest real X-ray 
source in the CDFs. These model PSFs were produced using 
the IDL-based s ource extrac tion tool ACIS_EXTRACT (for 
details see Bro os et alJl2004h . whereby the PSF for each in- 
dividual Chandra observation was calculated with the CIAO 
tool MKPSF using the CALDB PSF Ubraries, weighted by the 
number of counts in each exposure, and co-added. The near- 
est real source was always < 1' from the simulated source 
position. 
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10" 10' 102 iQj 10= 10" 10' 102 iQj 10-' 10= 

Measured 0.5-2,0 keV Counts Measured 2-8 keV Counts 

Fig. 2. — The simulated completeness (top) and flux recovery correction (bottom) functions for three representative positions in the CDF-N field as a function 
of measured counts. The solid, dotted, and dashed curves denote median values at off-axis angles of ()', 5', and 10', respectively, using all simulated sources 
within 2' of each representative position; the thin solid line indicates the value of unity. The soft- and hard-band corrections are shown on the left and right, 
respectively. Eddington bias is only significant for the 0.5-2.0 keV number counts in the 0' and 5' samples, where the recovered flux turns downward for low 
count sources. 




10" 10' 102 IQJ IQI 1Q5 IQO IQI 1Q2 IQJ ^ Ql ^ q5 

Measured 2-8 keV Counts Measured 2-8 keV Counts 

Fig. 3. — The individual flux recovery corrections for hard-band simulated sources within 2' of the 0' off-axis position shown in Figure|2](these sources are 
represented by the solid curve in the lower right diagram of Figure|2|. The solid curve indicates the median flux connection, akin to those plotted in the bottom of 
Figurel2l The corrections are shown both including (right) and excluding (left) Eddington bias, in order to demonstrate separately the effects of our photometry 
alone and our photometry plus Eddington bias. Photometry errors alone skew sources diagonally in these diagrams such that underestimated sources have both 
lower counts and higher flux recovery corrections, while overestimated sources have both higher counts and lower flux recovery corrections. Thus photometry 
errors by themselves imprint a gradual upward trend in the flux connections for decreasing source counts. The dashed curves indicate the sense of this photometry 
effect, showing tracks along which simulated sources with 10, 30, 100, 300. and 1000 counts would be scattered due to photometric enors. The tracks span a 
range of ±3(t. Very faint sources with low flux coirections tend to be simulated sources which lie below our nominal detection threshold but are detected because 
they sit on positive background fluctuations (both panels) or are caused by Eddington bias (right panel only). Adding Eddington bias to the simulation injects 
significantly more noise into the distribution of flux conections and provides several additional faint sources with low flux conections. As such. Eddington bias 
tends to pull the median flux coiTection down for faint sources. 
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Source searching and photometry of the simulated images 
were performed in a manner identical to that u sed to produce 
the CDF catalogs [i.e., using WAVDETECT jFreeman et al.l 
|2002) and custom software; see §3.2 and §3.4 of A03 for de- 
tails]. A completeness correction was determined by com- 
paring the number of simulated input sources to the num- 
ber of simulated detected sources as a function of detected 
counts. Likewise, a flux recovery correction was estimated 
by comparing the simulated input counts (both prior to and 
after including Eddington bias) to our measured aperture- 
corrected counts. Since both completeness and flux recovery 
vary across the CDF fields due to Chandra's changing PSF 
and spatially dependent vignetting, we determined the com- 
pleteness and flux recovery functions within a 2' radius region 
around each real source, averaged over the 200 simulations. 
The radius value of 2' was chosen as a compromise between 
the maximum area over which the PSF remains relatively 
constant and the minimum area needed to achieve reason- 
able statistics with our 200 simulations. The above correction 
functions were used to correct the measured source densities 
and fluxes, respectively. The completeness and flux recovery 
functions remain close to unity above 50-100 counts. Be- 
low this point, Chandra's varying PSF size and spatially de- 
pendent vignetting begin to affect source detection and pho- 
tometry. Figure 12] demonstrates how the completeness and 
flux recovery functions behave for three separate positions 
in the CDF-N field in both the soft and hard bands; similar 
functions were obtained for the CDF-S. Incompleteness is 
due to three factors: a few percent decrease for moderately 
bright sources due to occasional source overlap; a large de- 
crease for faint sources near the detection threshold; and a few 
percent increase for very faint sources due to source overlap 
or sources that lie in regions of particularly low background. 
The completeness curves shift as one moves off-axis because 
of the radially degrading sensitivity limit of the CDFs. Like- 
wise, deviations in the recovered flux correction are due to 
three factors: a few percent increase at all fluxes due to an 
additional aperture correction not originally accounted for in 
A03; a steady increase for faint sources due to photometry er- 
rors as detailed in Figure |3j and a steady decrease for faint 
sources due to Eddington bias. Eddington bias only appears 
to win out over photometry effects for soft-band sources near 
the centers of the fields. 

4. NUMBER COUNTS 

The cumulative flux distribution (logA^-logS') at each flux S, 
for all sources brighter than S weighted by the corresponding 
sky coverage, is 

N(>S)= Y^iCFi^iT'^ (1) 

i=Si>S 

where the sky coverage flj is the maximum solid angle over 
which each source could have been detected in both CDFs 
based on the 2<j limiting flux maps, and CFi is the complete- 
ness correction interpolated from the position- and count- 
dependent simulated completeness correction (see f|3|l. Addi- 
tionally, each flux S has been corrected for flux bias assuming 

Si = FRiS°, (2) 
where FRi is the position- and count-dependent simulated flux 
recovery function (see ^ and Sf is the original flux. Figure|3 
shows the total logA^-logS' in the soft and hard bands for the 
combined sample. The combined logA^-log5 curves are con- 
sistent with the distributions and models found by other au- 
thors (e.g., C02; M03) within the uncertainties. The CDFs 



appear to underestimate the soft-band number counts around 
^().5-2.o kev ~ 10"^"* erg cm"^ s~' compared to the model curve 
of M03, although this is only a « 2(t deviation. Figure|5lcom- 
pares the number counts in the CDF-N and CDF-S. We have 
worked out error bars on the cumulative distributions follow- 
ing Gehrels ( 1986) and calculated the deviation in quadrature 
at each data point in units of sigma. We find that the number 
counts from the two fields are consistent with each other at 
better than the la confidence level over the entire soft band 
and above F2-8 keV ~ 1 x 10"'^ erg cm"^ s"' in the hard band. 
Below this hard-band flux, statistical deviations gradually in- 
crease to 3.9(7 at the faintest flux levels, a finding similar to 
that presented in C02. This demonstrates that the field-to- 
field variations previously reported for the CDFs are entirely 
consistent with the lack of field-to-field variations found from 
ChaMP (Kim et al. 2004), as the latter study only examined 
the X-ray logA^-logS' for relatively bright sources where we 
find the CDFs are compatible. 

The total number counts in both bands have previously 
been best fitted with broken power laws with break fluxes 
between 10"'5-10""* erg cm"^ s"' (e.g., C02; M03). Given 
this break range and the limited statistics of the CDFs above 
~ lO"''* erg cm"^ s"', we estimate the slope only for the 
faint end of the number counts below 2 x 10"^^ erg cm"^ s"^ 
We have assumed a single power-law model of the form 
A^(> S) = A^i6(5'i/10~'*)~", where the slopes and normaliza- 
tions were determined with a maximum-likelihood algorithm 
(e.g. jMurdoch et alJfT973i) using the sky-coverage corrected 
(i.e., CFi^li), differential flux distribution. The best values of 
a were estimated by minimizing 

L = Mlna-M In ^(C/^i),)(5p - SJf{) 

i 

+ Y^n(CFini)-(a+l)Y^nSi, (3) 

where M is the number of sources used. The best values of 
Ni(, for a particular a were estimated by 



j:^(CFini)(sr-si_^o' 

The best-fit slopes and accompanying normalizations are pro- 
vided in Table[2 We typically do not have sufficient statistics 
to ascertain whether a single power law appropriately repre- 
sents the underlying flux distribution; however, visual inspec- 
tion of Figure |4] and the analyses below suggest that such a 
simple model is unlikely to represent adequately sources near 
or below our detection threshold due to the varied nature of 
the contributing source populations. 

Our best-fit slopes to the total faint-end number counts are 
0.55!°;!]^ (soft) and0.56!jj (hard). These values are slightly 
lower than previous estimates in both the soft (e. g., 0.70ib0.20 
from Mush otzkv et a l. 2000, 0.67±0.14 from |Brandt et alj 
2001b, 0.60±0.10 from Rosati et al. 200 i and 0.6 0!}}°? from 
M03) and hard (e.g., 0.61±0.10 from"^Rosati e t ah .2003 
0.63 ±0.05 from C02, and 0.44!;; from M03) bands, al- 
though they are consistent within measurement errors. Our 
lower values may be due to a number of the factors: (1) we 
are using data which probe a factor of « 2 deeper (i.e., the 
2 Ms CDF-N) where the number counts may become flat- 
ter, (2) we estimate the faint-end slope only using sources 
below 2 X 10"'^ erg cm"^ s"' (i.e., a factor of several lower 
than the break flux adopted by other studies), and (3) we use a 
somewhat different technique to calculate the incompleteness 
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Fig. 4. — The soft {left) and hai'd (right) number counts in the CDFs are shown in the upper panels for all sources (solid black curves) and three different 
X-ray subsets: AGN (solid and dashed grey curves), galaxies (solid and dashed dark grey curves), and stars (thick solid black curves). Also shown are the best-fit 
X-ray number count models of M03 (thin sohd black curves), the 1 Ms CDF-N fluctuation analysis results from Mivaii & Griffiths (2002, grey "fishtails"), 
and the pre di cted X-ray counts of star-forming galaxies from Ranalli et al. i2003), calculated ass uming the X-ray/radio correlation (e.g., Bauer et al. 2002^ 
|Ranalli et al. 2003) and the total (predominantly star-forming galaxies) radio number counts of Richards i2000, thin dotted grey lines). The percentage of the 
XRB resolved at a given X-ray flux (calculated from the number count models of M03, although renormahzed in the hard band to the total hard XRB flux density 
of lPe Luca & MolendiM2004h is shown at the top of the panels. Lower panels show the differential fraction of the total CDF soft- and hard-band samples each 
source type comprises. 



and Eddington bias corrections. The best-fitted number-count 
models from M03 are also shown in Figure|4]and display ex- 
cellent agreement apart from a slight soft-band deficit in the 
CDF number counts around f().5-2.o kev ~ 10"'"^ erg cm"^ s"' 
(presumably due to the effects of large-scale structure). At 
the CDF flux limits, we find total source densities of 9014;|;334 
and 5303^237 sources deg~^ in the soft and hard bands, respec- 
tively (see Table|2ji. 

4. 1 . Number Counts by Type 

A large fraction of the X-ray-bright CDF sources can be 
securely identified as AGN based on their X-ray luminosi- 
ties. X-ray spectral properties, radio properties, optical spec- 
troscopic classifications, and X-ray-to-optic al flux ratios (e.g. 
IXlexander et al. 2002at iBauer et alJ l2002al) . However, the 
classifications of many faint X-ray sources remain ambigu- 
ous because they lack firm redshifts, and their X-ray proper- 
ties are consistent with emission from either a low-luminosity 
active nucleus, star formation, or a combination of both. To 
understand these sources, we must make a few reasonable as- 
sumptions about the nature of the CDF sources as outlined 
below. Following these assumptions, we adopt two classifica- 
tion schemes, one which conservatively estimates the number 
AGN and one which conservatively estimates the number of 
star-forming X-ray sources. 



Optical magnitudes for the CDF sources have been mea- 
sured using Subaru observations in the CDF-N (Baraeret al] 
2003) and Wide-Field Imager (WFI) observati ons as part of 
an ES O deep public survey in the CDF-S (e.g.. lArnouts et all 
1200 Ih . Fifty-two of the CDF sources have no optical coun- 
terpart to the depths of the Subaru and WFI images {R > 27). 
We have assigned /?-band magnitudes to these blank X-ray 
sources assuming they represent the tail of the currently ob- 
served 7?-band counterparts. This assumption is validated by 
deeper HST imaging for a large subset of the CDF sources 
(e.g., Bauer et al. 2004, in preparation; iKoekemoer et alJ 
120041) . although we caution that a small fraction of sources 
could be much fainter than our estimates. Redshifts for 
the X-ray sources were culled from several recent sp ectro- 
scopi c (Baraer et al. 2003; Steidel et al. 2003; Le Fevr e et alJ 
'2004^'Szokolv et al. 2004; Wirth et al. 2004) and photometric 
( Alexander etal. 2001; Bai'seretal. 2003; Wolfetal. 2004 
catalogs. Of the 829 CDF sources, 425 have spectroscopic 
redshifts and another 80 have firm photometric redshifts. The 
remaining 324 sources that lack redshifts are nearly all opti- 
cally faint {R > 24) and have proven difficult to identify in 
large numbers even on 8-lOm class telescopes. In order to 
classify all sources, we first must estimate the redshifts for 
these remaining sources (thereby allowing X-ray luminosity 
and rest-frame absorption column determinations). We note 
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Table 1. Number Count Slopes and Normalizations 



(1) (2) 
Class Subclass 


(3) (4) (5) 
# Slope (a) Ni(, 
S H S H S H 


All Total 


724 520 


55+!!!!3 56i!-!1 


3039+?L 7403+'25 

9? 5')') 


AGN Total (opt.) 

Total (pess.) 

Optical Type 1 

Optical "Not Type 1" 

X-ray Unabs. 

X-ray Abs. 
Galaxies Total (opt.) 

Total (pess.) 

Starbursts (pess.) 

Quiescent Gal. 

Ellipticals (pess.) 
Stars Total 


599 484 
545 466 
73 65 
526 419 
Tin 1/10 

368 342 
157 30 
103 12 
55 5 
35 5 
13 2 
22 3 


o:47| o:48| 

A 47+0.04 37+""' 

"■^'-0.03 ""^ -0.10 
n ^n+O.OS n qc+0.20 

62+" "'' 44™'" 
"■°^-0.04 "-^^-0.11 

1 26+<'-"' 2 02+" " 

-fl-Ofi -».3I 

1 oo+fl.IO 9 CI +0.53 

1.70+^!= — 

-0.12 

^■^^-0.31 

6C+0.12 _ 


2625+'^^ 6901!^'* 

7365+'** 6749+'" 
-104 °'^'-201 

427+150 

^^'-103 

2266!«« 4207^495 

QQ8+3^9 1 Q^O+444 

"°-15T -314 
1578+g 4337+{{4 

438!" 4081!i358 

263^1 3673igl 

108+? — 

-12 

115+'3 
56+i^ 
125+11 


log(Lx) > 44.5 

43.5^4.5 
42.5^3.5 
41.5^2.5 
40.5^1.5 


23 24 
165 169 
276 229 
115 57 

17 4 


1 1+" "2 44+0-2S 

43+0.04 31+011 

-Q.Q5 -Q.13 

1 29+0.05 , 45+0.20 
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-0 ''3 


780+^"^ 1553+2QQ 
1380+1"* 2533+]28 
371+^^ 4026+1'*' 
4S+1 — 


log(iVH) 23-24 
22-23 
21-22 
<21 


100 142 
250 180 
145 103 
724 55 


OQft+O.BIi 33+^'** 
"■'':'-0.08 "■-'-'-0.33 

51+" "5 60+<''* 
-QQ5 -Q.12 
10+0 05 98+<'-2- 

Q 54+0.69 , ] 4^0.13 
"■■^ -0.08 -0.35 


38i:j^ 1619«^ 
1097^^^ 2974+3« 
5f:4+204 2231+''* 
-'°^-163 " -411 

483!^^ 985!=™ 



Note. — T he two entri es in c olumns 3-5 denote the soft (S) and hard (H) band results. Column 1: Source classification (see i|4.1.U . Column 2: Source 
subsample (see i|4. 1.2l and i|4. 1.41 . Column 3: Number of sources in class in the soft and hard samples, respectively. Columns 4 and 5: Number-count slopes 
(a) and normalizations (Nk,, sources deg - at lO""" erg cm ^ s ') as estimated using the maximum likelihood method (e.g.. lMurdoch et ali, 1973.1 on the soft 
and hard samples, respectively. Fitting was performed using only the sources with fluxes below 2 X 10^" erg cm~^ s~' (i.e., below the known break in the X-ray 
number counts) for all classes except that of "Stars". For some cases, there were too few sources or too much scatter below 2 X 10"'^ erg cm"^ s"' to provide 
reliable number-count parameters. Such cases are denoted by " — ". 
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Fig. 5. — Comparison of CDF-N (solid curves) and CDF-S number counts 
(dashed curves) for the soft (lower set) and hard (upper set) bands. The two 
fields are consistent at the Icr confidence level over the entire flux range sam- 
pled in the soft band and above F2-8 keV ~ 1 x 10~'^ erg cm~^ s~' in the hard 
band. 



that detailed studies of optically faint X-ray sources indicate 
that they often have hard X-ray colors typical of obscured 
AGN and red optical colors typical of early-type galax ies at 
z ^ 1-3 (e.g., Alexa nder et al.l200lH Ko ekemoer et al.ll2004t 
iTreister et al. 2004). 

In Figure |6| we show the distribution of redshift versus R- 
band magnitude for the CDF sources with known redshifts, 
as well as templates for an Mb = -23 QSO and unevolved 
M* Sc, Sb, and E galaxies. /iT-corrected galaxy tracks were 
derived from the study of Poggianti (1997). Since no evo- 
lutionary corrections have been made, these tracks should be 
considered extreme in the sense that a typical z = 1 ellipti- 
cal should have bluer colors than those shown here. The K- 
corrected QSO track was calculated using a custom compos- 
ite QS O spectrum (consistent with that of Vanden Berk et ak] 
I2OOIL where they overlap) assuming the QSO continuum has 
a = 0.5 (where /v v'"'), typical emission-line strengths, and 
standard absorption due to the Lyman alpha forest. 

While the broad-line AGN span a wide range in R and are 
generally consistent with the QSO track, the overwhelming 
majority of the other CDF sources follow the galaxy tracks, 
indicating that their optical light is likely to be dominated by 
their host galaxies (see Grogin et al. 2003 for confirmation of 
this dominance of the optical host galaxies from HST imag- 
ing). We note that the correlation found by Fiore et al. (2003|) 
between intrinsic Lx and fx/ fo can essentially be reduced to 
one between optical magnitude and redshift (such as is seen 
in Figure |5) as long as we include a correction for X-ray ab- 
sorption, which these authors suggest is typically <15% for 
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(1) 

Class 



All 



AGN 



Galaxies 



Stars 



log(Lx) 



(2) 
Subclass 



(3) 

Total # deg" 



Total 



Total (opt.) 
Total (pess.) 
Optical Type 1 
Optical "Not Type 1'' 
X-ray Unabs. 
X-ray Abs. 
Total (opt.) 
Total (pess.) 
Starbursts (pess.) 
Quiescent Gal. 
Ellipticals (pess.) 
Total 



>44.5 

43.5-44.5 

42.5^3.5 

41.5^2.5 

40.5-41.5 



log(A'H) 23-24 
22-23 
21-22 
<21 



(4) 

CDF XRB % 
S H 



(5) 

TOTAL XRB % 
S H 



9014 



+347 
334 



7166!^'" 
6342 

325«; 
6840+311 
1774!™ 
5392!!' 

2552 



27| 
38 



\121_ 
^*^-132 



343 
399: 



+1^ 



4558!- 
3970+11 
350!*f 
4207!|5 

1430!™ 
3J27+178 

1298+11 
7 11+™ 
123!| 

481+327 
-202 
105+'3' 

34+f3 



109^ 

868^3 

2716+113 



3055: 
411 



"3! 

1050i^ 
2157+152 

1168+!i 



1754+|« 

862!:*] 
1098!|i« 



58 



- --28 

1668+!" 
1199!]30 
435!«2 



70.3!^;^ 



64.6^; 
33.8!t^ 
32 2+2'' 
47.8!^i 
18.1!';3 
q 9+0.3 

qqJ:? 

"■^-0.1 

o.6+<;' 



75 4+5-3 



':^:i 



0.2 
2.5 



13.0^ 
33.7!|^ 
15.9!'3 
3.1 
2.6 



73.0: 
71.5 

20.6!|'j 

52.4!3:8 

29.8!3] 

43.2!3:3 

2 6^-^ 
-Q-5 
1 2^<).4 

4+^'^ 
"■^-0.2 

7+" '* 

-0.3 

0.1 
1.3 



^:i 



1+U.4 
-"-0.3 

14.3+1? 



12 7+-+2 
33 9+3 3 
22,1+1.9 

4 2+^:^ 

^■^-0.6 
I 7+1.4 
^■'-0.8 



3.0+ 



17.2+ 



1 

30.2 



4+1.9 
^-1.7 
+3.6 



-1.7 

23.9!2-3 

14.8!|| 
16.7!" 



89.5!;^;^ 



86.9!" 



74.6^ 
73.2+5-9 

41 5+^-^ 
33.1+'*-5 
55.5+1^ 
19 l+'* 2 



83.2 
81 

28.5 



37 



4.1+ 



i+4;o 

-4.0 

2.7 (< 6.7) 



I n 2+4. 1 
"'■■^-4.0 



^.8 
54.7!5-5 

7+5:1 

45.5!'-? 
2.8 (< 6.8) 
1.4 (< 5.6) 



1.4 (< 5.6) 



Note. — T he two entr ies in columns 3-6 denote the soft (S) and hard (H) band results. Column 1: Source classification (see i|4.1.11 . Column 2: Source 
subsample (see i|4. 1.2l and > I4.1.4I . Column 3: Source number density in the soft and hard samples, respectively. Column 4: Percentage of the XRB flux density 
which each source class contributes (considering detected CDF sources only). Column 5: Percentage of the XRB flux density which each source class contributes 
(considering CDF sources and sources brighter than those in the CDFs, which we estimate to comprise a total of 19.2% in the soft band and 11.5% in the hai'd 
band; see ||4]. For some classes, the percentage of the XRB for brighter sources is not well known and has been omitted (denoted by " — "). 



w90% of their so urces and rarely ex cedes factors of a few. In 
this new form, the lFiore et alJ i2003l) correlation tracks bright 
optical sources as well as the galaxy templates, but appears to 
overestimate redshifts for sources with 7? > 24 by increasingly 
large factors. Thus caution should be exercised when extrap- 
ola ting this correlatio n to extremely faint optical magnitudes 
(cf. lPadovani et alJ20 04). 

The fainter sources that lack redshifts are likely to con- 
tinue to follow the galaxy tracks in Figure |5J and thus in 
the absence of reliable photometric redshifts we can use this 
good correspondence as a crude redshift indicator to glean ba- 
sic trends. Differences between the galaxy templates do not 
appear large below R ^ 27. However, the Sb template for- 
mally provi des the best empirical fit to the available data (see 
also, e.g., Alexander et al. 2002b; Barger et al. 2002', 2003) 
and was therefore adopted as our median spectral energy dis- 
tribution (SED) to convert /?-band magnitude to redshift. We 
assumed a scatter about this template equal to that measured 
between R=23-24 for the known redshift distribution. We 
caution that this technique is only valid in a statistical sense 
and is only strictly true if (1) fainter X-ray sources remain 
host-galaxy dominated at optical wavelengths and (2) their 
typical hosts are well-represented by an Sb galaxy SED. Re- 
assuringly, the small number of optically faint sources with 
redshift determinations follow this template. 

We must also calculate rest-frame X-ray absorption column 
densities for all CDF sources in order to determine the intrin- 
sic power of the AGN. These were determined from direct X- 
ray spectral analysis (Bauer et al. 2004, in preparation), which 
is more accurate than using simple band ratios. An absorbed 



power-law model (wabs+ ;?o: lMorrison & McCammoiJl983h 
was fitted to all CDF sources. The absorption and normaliza- 
tion were varied as free parameters. The photon index of the 
power law was allowed to vary such that F > 1.7 for sources 
with more than 100 counts in the 0.5-8.0 keV band, while 
it was fixed at F = 1.7 for sources below 100 counts. The 
fitting was perf ormed on the unbinned X-ray data using the 
Cash statistic (iCashlll97^ to maximize spectral information 
for low-count sources. The spectral parameters of the CDF 
sources were then used to calculate unabsorbed, rest-frame 
0.5-8.0 keV luminosities. 

4.1.1. Source Classification 

Our first source classification scheme is based on intrin- 
sic 0.5-8.0 keV luminosities. X-ray spectral properties, radio 
properties (slopes, morphologies, and variability), and opti- 
cal spectroscopic classifications. Our specific AGN criteria 
are motivated by the following findings. From X-ray obser- 
vations in the local Universe, purely star-forming galaxies do 
not appear to have X-ray luminosities exceeding Lq, 5-8.0 keV ~ 
3 X lO'*^ erg s~', and rarely, if ever, have intrinsic absorp- 
tion column densities abo ve A^h ~ 10^^ cm"^ (e.g., Fabbian^ 
1989; Colb ert et a D|2004). The above numbers are proba- 
bly conservative, considering the most X-ray-luminous local 
star-forming galaxy known, NG C 3256, only has an X-ray lu- 
minosity of Lo.5-8.0 keV < 10'*2 (iMoran et alJll99^ iLira et all 
2002) and that X-ray emission from star-forming galaxies is 
typically extended; obscuring X-ray emission with an average 
absorption column density of A^h ^ 10^^ cm"^ would require 
extraordinary amounts of intervening gas. Finally, nearby 
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Fig. 6. — The R-band magnitude versus redshift distribution for all 829 CDF sources, divided into sources w itli known redshifts (left) an d unidentified CDF 
sources (right) adopting an R-band-to-redsliift estimate based on an Sb galaxy spectral energy distribution (see i|4. 1 . 1 1 'Colem an et a 111980) . The sources with 
known redshifts have additionally been subdivided into spectroscopically (Z) and photometrically (PHOT Z) identified objects, as well as those with broad 
emission lines (BLAGN). Templates for an Mg = -23 QSO and unevolved M* Sc, Sb, and E galaxies are shown overlaid, adopting A'-corrections based on a 
custom composite QSO spectrum (consistent with that of Vanden Berk et al. 2001 , where they overlap) and the galaxy models of Poaaianti ( 1997), respectively. 
For comparison, we also plot the correlation between intrinsic Lx and fx/ fo discussed by Fiore et al. ( 2003), which can be reduced to a correlation between 
optical magnitude and redshift provided we apply a !5il5% correction factor to account for the typical intrinsic X-ray absorption these authors find in Lx (for 
details see .Fiore et al, 20 03V Additionally, we show corresponding R-band and redshift histograms to illustrate the overall distributions. The dotted histogram 
in the R-band magnitude distribution represent the 52 sources with no known redshift and no R-band counterpart (see i|4. ll for details). The dashed curve in the 
redshift histogram indicates the expected distribution of AGN measured from R05/ir observations iGilli 2003), normalized to the total number of AGN in the 
CDFs (assumed to be ~80% of all CDF sources; see column 3 of Tablel2l. The large discrepancy between ROSAT and CDF sources at low redshift is due to the 
fact that CDF AGN appear to peak in number density at somewhat lower_ redshifts (e.g., Ba raer e t al-200.3;.Ueda et al..200.3il . There is also a non-negligible and 
increasing fraction of star-forming galaxies in the CDF sample (e.g., iHornschemeier et ali20Q3l) . 



extragalactic X-ray sources with extremely flat X-ray spec- 
tral slopes (r < 1) are almost always identified as highly ob- 
scured or Compton-thick AGN, where the primary emission 
from the AGN is almost completely obscured and we only see 
the flat scattered or reflected component (e.g..'Ma iolino et aH 
[1998; Bassani et al. 1999). At radio wavelengths, the vast 
majority of sources with powerful radio jets/lobes, flat ra- 
dio spectral slopes, or strong radio variability are AGN (e.g., 
[Condon 1984, 1986), while extragalactic sources with broad 
(EW>1000 km s~') or high-excitation optical emission lines 
are almost universally classified as AGN (e.g., Osterbrock 
[l989). 

Based on the above constraints, the CDF sources were di- 
vided into AGN, galaxies, and stars. Twenty-two sources 
were classified as Galactic stars based on their spectroscopic 
identifications. Six hundred thirty-two sources were classi- 
fied as AGN based on at least one of the following properties: 
Na > 10^^ cm"^, hardness ratio > 0.8 (equivalent to effec- 
tive r < 1.0), Lo.5-8.0 kev > 3 X 10"^^ erg s"', or broad/high- 
ionization AGN emission lines. We recognize that there 
will be AGN-dominated sources not selected by these crite- 
ria, such as AGN with F > 1.0 and rest-frame luminosities 
io.5-8.0 keV < 3 X lO'*^ erg s"' ; however, such AGN are difficult 
to classify even locally. Seventy-six sources were classified as 
star-forming galaxies based on either having off-nuclear X- 
ray emission or A^e < 10^^ cm"^, hardness ratio < 0.8, and 
Lo.s-a.o keV < 3 X 10'*^ erg s"' . The classifications of a remain- 
ing 99 CDF sources were considered ambiguous due to weak 
X-ray spectral constraints; we have tentatively classified these 



sources as star-forming galaxies, although there is likely to be 
some degree of obscured and low-luminosity AGN contami- 
nation. Note that the numbers provided above are based on the 
classification of all 829 CDF sources, and should not be con- 
fused with the numbers given in Tabled which indicate the 
number of each source type used in the soft-band and hard- 
band number-count estimates, respectively. 

Figure compares 7? -band magnitudes and 0.5-8.0 keV 
fluxes for all 829 CDF sources, separated into several un- 
absorbed, rest-frame 0.5-8.0 keV luminosity classes. Those 
sources classified as AGN under our first scheme are high- 
lighted. From inspection, it is apparent that the X-ray-to- 
optical flux ratio crudely tracks X-ray luminosity (and hence 
AGN activity), explaining why fns-s ji kpy/fw =0.1 is a use- 
ful AGN/galaxy discriminator (e.g., Maccacaro et al. '198^ 
Stocke et al. 1991; Schmidt et al. 1998; Akivama et al. 200(1 
IHornschemeier et ab 2000). At very faint X-ray flux lev- 
els, however, the relation will begin to break down for cer- 
tain source types due to the very different nature of their 
X-ray and optical K-corrections (most notable in this regard 
are Type 2 Seyfert galaxies, whose X-ray-to-optical flux ra- 
tios can vary by more than an order of magnitude between 
z = 0-2; Moran et al., in preparation). There is also likely to 
be some degree of obscured and low-luminosity AGN con- 
tamination at /0.5-8.0 kev/./R < 0.1, as well as considerable 
potential galaxy contamination for /0.5-8.0 kev//R > 0.1 and 
■P{).5-8.o keV ^ 4 X 10"'^ erg cm"^ s"'. We attribute the fact 
that star-forming galaxies begin to have /0.5-8.0 kev//R > 0.1 
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Fig. 7.— R-band magnitude versus 0.5-8.0 keV flux for all 829 CDF 
sources. Sources are divided into several estimated, unabsorbed 0.5-8.0 keV 
luminosity bins (color-coded following the order of the rainbow, with red 
denoting the highest luminosity and a bsorpt ion bins). Sources that are con- 
sidered AGN following the criteria in i|4. l.ll are denoted by small white dots. 



at faint X-ray fluxes to (1) a stronger K-correction in the op- 
tical band than in the X-ray band that tends to shift potential 
star-forming galaxies into the AGN region and (2) increas- 
ingly large uncertainties in the source band ratios (the primary 
AGN discriminator for low-luminosity X-ray sources) at faint 
fluxes. 

We can estimate the prevalence of obscured AGN among 
these ambiguous faint X-ray sources by stacking them in suit- 
able subsets and examining their average X-ray band ratios. 
Given the good correspondence between X-ray-to-optical flux 
ratio and AGN activity, we stacked all of the ambiguously 
classified X-ray sources with unconstrained band ratios and 
■Pb.5-8.0kev < 4 X 10""' erg cm~^ s~' in decades of decreas- 
ing X-ray-to-optical flux ratio. This yielded average effec- 
tive F's of 1.35, 1.45, 1.77, and 2.01 (associated errors are 
±10-20%) for /0.5-8.0 kev//R = 1-10, 0.1-1, 0.01-0.1, and 
< 0.01, respectively. The steep effective photon indices for 
sources stacked in the "galaxies" region of Figure sug- 
gest that there are relatively few obscured AGN among these 
sources and that they are likely star-forming in nature (al- 
though we cannot exclude the presence of soft low-luminosity 
AGN or AGN/starburst composites). The flatter effective pho- 
ton indices for sources in the traditional AGN region, on the 
other hand, suggest that a significant fraction are in fact ob- 
scured AGN, with the rest presumably a mix of unobscured 
AGN and star-forming galaxies. Thus, our first classification 
scheme appears to provide a pessimistic estimate of AGN and 
an optimistic estimate of star-forming galaxies (since some of 
the ambiguous sources classified as star-forming galaxies may 
in fact be powered by AGN). 

Our second classification scheme is a sli ght variation of 
the first and is similar to that presented in Alexand er et al.l 
( 12002 a) and .Bauer et al, (2002aJ . In addition to the X-ray 
spectral properties, intrinsic X-ray luminosities, radio mor- 



phologies, variability, and optical spectroscopic classifica- 
tions, we further classify ambiguous sources based on their 
X-ray-to-optical flux ratios. Thus we considered 698 sources 
to be AGN based on at least one of the following proper- 
ties: fo..5-8.(i kev / fR > 0.1, A^H > 10^^ cm~^, hardncss ratio 

< 1.0), L0.5-8.0keV > 3 X 



/0.5-8.O kev//R ■> U.i, JVh 

> 0.8 (equivalent to effective F 

-l42 e-1 



10 erg s" , or broad/high-ionization AGN emission lines. 
The remaining 109 sources were considered star-forming 
galaxies, although again some low-luminosity AGN might be 
present. In contrast to our first scheme, this division should 
provide a more optimistic estimate of AGN and a more pes- 
simistic estimate of star-forming galaxies, as it attempts to ad- 
dress the issue of strong AGN contamination for unclassified 
/0.5-8.0 kev//R > 0.1 sources. 

Number counts for the two classification schemes are com- 
pared in Figure|3 and our best-fit slopes to the total faint-end 
number counts are provided in Table^for each category. Im- 
portantly, the slopes of the AGN and galaxy number counts 
are not strongly affected by the adopted scheme. The true 
number counts for each class are likely to lie somewhere in 
between these two determinations. Given that the stacked ef- 
fective photon indices of the ambiguous sources in the AGN 
region are relatively hard, and that there is likely to be some 
soft AGN contamination in the "galaxies" region, we will 
adopt the second scheme presented above (optimistic AGN, 
pessimistic galaxies). 

We can compare our pessimistic galaxy sample to the 
Bayesian selected sample of normal/starburst galaxies in the 
CDFs from Norman et al. (2004). Using only CDF sources 
with spectroscopic redshifts z < 1.2 and a selection based 
on X-ray luminosities, observed X -ray har dness ratios, and 
X-ray-to-optical flux ratios, iNorman et alJ (j2004) identified 
a total of 210 normal/starburst galaxies. This number is 
substantially larger than the 109 galaxies in our pessimistic 
galaxy sample or even the 175 galaxies in our optimistic 
sample and warrants investigation. Part of this discrepancy 
arises from the fact t hat Nor man et alJ (2004) use the cata- 
log of iGiacconietalj (200^ rather than that of A03; there 
are 18 iGiacconi et alJ ( l2002h sources in their sample which 
are not in the main catalogs of A03. Of the remaining 192 
Norman et al. (2004) galaxies, we find 120 and 94 matches to 
our optimistic and pessimistic samples, respectiv ely. Of the 
55 sources classified as galaxies here but not by N orman et alJ 
(2004), ^ 30% were objects lacking spectroscopic redshifts 
or with z > 1.2 (15 cases) while it is unclear why the rest 
failed their criteri a (40 cases). Conver sely, sources classi- 
fied as galaxies by 'Norman et al. (2 004h but AGN here were 
objects which had a combination of /0.5-8.0 kev//R > 0.1 (62 



cases), A^H > 10 cm (44 cases), hardness ratio > 0.8 (23 
cases), or Lo.5-8.0 keV > 3 x 10"*^ erg s"' (20 cases). Based on 
our criteria, we exp ect sign ificant AGN contamination to be 
present in the lNorman et al.l ( l2004l) sample. 

We note that the upward trend seen below Fo.5-2.0 keV 5 x 
10"'^ erg cm"^ s"' in the AGN number counts in Figure |4] 
and elsewhere is likely caused by obscured AGN which have 
composite X-ray spectra, such that a significant fraction of the 
soft-band flux is likely due to either star formation associated 
with the host galaxy or complex AGN spectra (e.g., partial 
covering, scattered radiation, or reflection). 

The AGN number-count slopes are flat in both bands, indi- 
cating that we are now observing AGN beyond the differential 
peak in the AGN number-count distribution. While the X-ray 
source densities of AGNs, in a differential sense (see lower 
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panel of Figure appear to be in a state of decline at X-ray 
fluxes below ~ 10"'^ erg cm"^ s"', the star-forming galaxy 
population is rising strongly (nearly Euclidean and consis- 
tent with the previou sly measured norm al galaxy X-ray num- 
ber counts of Horns chemeier et aD '2003). The star-forming 
galaxy number counts also a gree well with the predictions 
made by Rana lireit alJ J200 3I) using the X-ray/ra dio correla- 
tion (e.g., Bauer et al.'*2002aVRanalli et al.l l2003h and the ra- 
dio number counts of Richards (2000), indicating that there is 
a good correspondence between the X-ray and radio emission 
in these sources. 

For completeness we also present the X-ray number counts 
for the 22 spectroscopically identified. Galactic stars discov- 
ered in the CDFs. These high-latitude stars are typically old 
G, K, and M types which are thought to emit X-ray emission 
via magnetic flaring; we refer the reader to Feigelson et al. 
(|2004), who have performed analyses on a subset of the X- 
ray-detected stars presented here. The number-count slope of 
the Galactic stars is flat in the soft band and indeterminate in 
the hard band. 

4.1.2. AGN Number Counts 

We further subdivided the AGN sample based on both opti- 
cal and X-ray properties. Unfortunately, the only meaningful 
optical division we can place is based on AGN with broad 
lines reported in their optical spectra (i.e.. Type 1), and those 
apparently lacking them (i.e., not obviously Type 1 AGN), 
since the optical spectra are not archived, adequately ana- 
lyzed, nor typically of very high signal-to-noise. Thus the 
"not obviously Type 1" category is probably too loose to pro- 
vide much insight, as it is likely to include not only clas- 
sic narrow-line AGN (Type 2), but also slightly more lumi- 
nous versions of Compton-thick AGN like NGC 4945 and 
NGC 6240 (e.g.. Matt et al. 2000), XBONGS JComastri et alJ 
|2003), and other odd types. We further caution that spec- 
tral misclassifications are likely for some sources since broad 
lines may be apparent in certain parts of the optical spec- 
trum but not others or may be missed due to host-galaxy con- 
ta raination or POor si gnal-to-noise data (e.g., see discussion 
in Moran et al. 2002). In particular, there are a significant 
number of CDF sources with high photometric redshifts that 
could turn out to be broad-line AGN (see Figure |6j. Given 
this likely incompleteness, the number counts for the optical 
Type 1 and not obviously Type 1 AGN must be considered 
lower and upper limits, respectively. Similarly, AGN with 
column density estimates A^h < 10"^ cm"^ were considered X- 
ray unobscured/mildly obscured AGN, and those above X-ray 
obscured AGN. Again, we caution that this distinction is far 
from exact, as it is based on the fit of a relatively simple spec- 
tral model, which surely is inadequate considering the known 
spectral complexity exhibited by local AGN (e.g., reflection, 
scattering, partial covering, host-galaxy contamination). Such 
effects tend to contribute more soft flux than predicted by a 
simple model and therefore can lead to an underestimate of 
the column density or confuse the spectral fit. 

From Figure |8] we find that optical Type 1 AGN fall off 
abruptly below ~ 10"'"^ erg cm"^ s"' in both bands, ulti- 
mately reaching source densities w 10-20 times lower than 
those of AGN that are not obviously Type 1 at the CDF de- 
tection threshold. While this may be a legitimate effect, the 
ability to detect broad lines tends to decrease at fainter op- 
tical (and hence X-ray) fluxes, and thus we cannot exclude 
the possibility that observational limitations contribute to this 
decline. Considering AGN by their X-ray properties yields 



similar, albeit less extreme, results. The source density of 
unobscured/mildly obscured AGN gradually trails off below 
^ lO"''* erg cm"^ s"' in both bands, ultimately reaching source 
densities ss 2-3 times lower than those of obscured AGN at 
the CDF detection threshold. 

The faint-end number-count slopes of all AGN subsamples 
are < 1 . At the CDF flux limits, we find AGN source densities 
of 1\66^^\ sources deg"^ (soft) and 4558^20? sources deg"^ 
(hard), which are factors of ^ 10-20 higher than the 300-500 
sources deg"^ found in the deepest optical spectroscopic AGN 
surveys (e.g., Hafl et al. 2000; Steidel et al. 2002; Wolf et alj 
2003; Hunt et al. 2004). At face value, this means that optical 
observations appear to be missing >90% of AGN compared 
to the deepest X-ray observations, although we concede that 
a direct comparison is not entirely fair. The optical surveys 
often target specific redshift and luminosity ranges, and X- 
ray and optical survey results are more consiste nt when com- 
pared within these ranges (c.f. jHunt et alJ20'04l) . Importantly 
though, the AGN discovered in deep optical surveys are gen- 
erally AGN with strong emission lines or blue continua and 
only have source densities comparable to our optical Type 1 
AGN sample alone (i.e., «325-350 sources deg"^). Given 
these facts. X-ray observations therefore appear to be much 
more effective at identifying AGNs, particularly (more typi- 
cal) obscured AGNs. This conclusion is not entirely surpris- 
ing given that the optical spectra of many AGN in the CDFs 
are dominated by host-galaxy light. 

4.1.3. AGN Completeness 

While the AGN source densities achieved above are im- 
pressive, it is important to understand the completeness of 
our X-ray selection of AGN. It is likely, for instance, that 
we are still missing a significant fraction of low-luminosity 
AGN (primarily those which have steep X-ray spectra or fail 
to dominate over host-galaxy emission), as such AGN are 
often difficult to identify even in the local Universe (e.g., 
iHo etaUfTool . Local studies have also shown that ^40% of 
moderate-to-high luminosity AGN are Compton thick, with 
very little direct emission from the AGN (e.g., Maiolino et afl 
1998; Mattetal. 2000). Again, such sources are difficult 
to identify since their observed X-ray luminosities are often 
i'0. 5-8.0 kev < 10"*^ erg s"' and contaminated by host galaxy 
emission; many such sources may be present in the CDFs but 
fall below our AGN luminosity threshold. 

Putting these caveats aside, however, there are only two 
classified AGN in the CDFs which are not yet detected 
at X-ray wavelengths (both in the 2 Ms CDF-N). The 
first, VLA J123725.7+621128, is a radio-bright (w 6 mJy 
at 1.4 GHz) wide-angle-tail source estimated to lie between 
z ~ 1-2 with a rest-fram e 0.5-2.0 keV luminosi ty limit of 
< 5 X lO"*' erg s'^ (e.g., 'Snellen & Best .2001; .Bauer et alJ 
2002b).'^ The sec ond, 12 3720.0+621222, is a narrow-line 
AGN at z = 2.445 f Hunt et al.ii2004 with a rest-frame 0.5- 
2.0 keV luminosity limit of < 2 x 10"*^ erg s"' . 

There are several more AGN candidates which lack X- 
ray detections, although their classifications are considerably 
more tentative. Most not ably, 10 of the 30 radio sources with 
-1-10 mJy at 1.4 GHz (IRichar ds"2000^. A. Koekemoer 2004, 
private communication) lack X-ray detections in the com- 
bined CDFs. This number is consistent with expectations that 
—30% of radio sources at 1 mJy are star-forming galaxies 

VLA J123725. 7+621 128 can be marginally detected with some tweaking 
of WAVDETECT's detection parameters beyond the scope of A03. 
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Fig. 8. — Plot identical to Fig.|4]except we plot here only the number counts for the total (solid black curves) and subsets of our AGN class of objects: 
total AGN (thick solid black curves), optical Type 1 AGN (dashed light grey curves), optical "not obviously Type 1" AGN (solid light grey curves), X-ray 
unobscured/mildly obscured AGN (dashed dark grey curves). X-ray obscured AGN (solid dark grey curves). 



(e.g., LTackson & Wal]lfT999h . although four of these sources 
have faint (R > 23) optical counterparts (atypical of compa- 
rable radio-detected star-forming galaxies) and thus may be 
radio-loud AGN. At optical wavelengths, Wolf et al. (2004) 
report 51 "QSO" candidates within the 1 Ms CDF-S region 
using the COMBO- 17 photometric redshift survey. Nearly 
half (24) of these lack X-ray counterparts. All but one of 
these non-detections, however, have ambiguous classifica- 
tions [i.e., "QSO (Gal?)"] or R > 23.5 (where the rehabil- 
ity of the COMBO-17 photometric redshifts becomes poor). 
Thus, we do not consider these to be convincing. Finally, 
ISaraiedini et al.l(l2003i) discovered 16 (with a correction of ±2 
for incompleteness and potential spurious sources) variable 
galactic nuclei within the Hubble Deep Field-North. These 
variable nuclei are thought to be primarily low-luminosity 
AGN (-15 < Mb < -17; e.g.. Ho et al. 1997), although a few 
may be nuclear supernovae (e.g., Riess et al. 1998). Only 
six have X-ray detections in the 2 Ms CDF-N, indicating 
<40% overlap with the X-ray sample. If all of the vari- 
able nuclei are indeed AGN, then their number density ('^ 
1 1000 ± 2000 deg-2) would surpass that found in the CDFs. 
It should be strongly emphasized here, however, that many of 
the AGN selected via optical nuclear variability are likely to 
have X-ray emission that is intrinsically 1-2 orders of mag- 

^ ISara iedini et al. 120031) actually report 7 X-ray matches, but the match of 
CXOHDFN J123651.73+621221.4 to HDF 3-461.9 is incon-ect as this X- 
ray source lies outside of the galaxy and has a more viable faint near-IR 
counterpart. Also, the X-ray luminosities for several of the matches are low 
enough to be due entirely to star-formation rather than AGN activity. 



nitude lower tha n the CDF dete ction limits (i.e., ^ 10^^- 
10^*^ erg s~'; see iHo etalJ 12001 . and thus are not directly 
comparable. 

Given all of the above. X-ray emission appears to be one 
of the most efficient and complete selection criteria currently 
available for selecting moderate-to-high luminosity AGN. 

4. 1 .4. Galaxy Number Counts 

As noted above, the star-forming galaxy population is ris- 
ing rapidly below 10"'^ erg cm"^ s"' in both bands. As 
was done for the AGN, we have further subdivided the galaxy 
sample based on optical and X-ray properties to gain fur- 
ther insight into this popul ation. We used the //5r-derived 
elliptical galaxy sample of llmmler et alJ J2004h to separate 
ellipticals (presumably passively evolving at the redshifts of 
relevance here) from spirals/irregulars/mergers (presumably 
actively star-forming), on the assumption that the former is 
likely to be dominated by different X-ray emission mecha- 
nisms (e.g., low-mass X-ray binaries, low-luminosity AGN) 
compared to the latter (e.g., high-mass X-ray binaries, su- 
pernovae, hot gas). The Immler et al. (2004) elliptical galaxy 
sample is derived from the Great Observatories Origins Deep 
Survey (GOODS) observations, which cover only w40% of 
the CDFs, so the ratio of ellipticals to other types is probably 
underestimated. Fortunately, 82% of the galaxy sample over- 
laps with the GOODS regions, so the elliptical galaxy num- 
ber counts are likely to be underestimated by < 10%. The 
non-ellipticals were split into two groups above and below 
/0.5-8.0 kev //r = 0.01, as set forth bv lAlexander et alJ J2002ah 
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Fig. 9. — Plot identical to Fig.l^except we show here only the number counts for the total (solid black curves) and subsets of our galaxy class of objects: all 
galaxies (solid dark grey curves), starburst galaxies (dashed grey curves), quiescent galaxies (solid light grey curves), and el liptical galaxies (thick solid black 
curves). Open symbols indicate the stacking results for (1) Lyman Break galaxies (Brandt et al. 2001a), (2) spiral galaxies IHomschemeier et all20O2l) . (3-4) 
Lyman and Balmer Break galaxies tNandra et al. 2002) . and (5) very red objects {I—K > 4: .,Alexander et al. 2002b). 



to denote starburst and quiescent galaxies, respectively. 

The number-count slopes of the starburst and quiescent 
galaxy samples are consistent with being Euclidean (see Ta- 
bleQ, suggesting they may track to some extent the strongly 
evolving infrare d galaxy population (e.g., Elbaz et al. 1999; 
[Alexander e t al. 2002a). The elliptical galaxies, meanwhile, 
have a somewhat flatter slope. Caution should be exercised 
when interpreting the very steep slopes observed in the hard 
band, of course, as these are uncertain and may be biased up- 
ward by unidentified low-luminosity or Compton-thick AGNs 
in all three samples. Quantifying the star-forming galaxy 
number counts at brighter X-ray fluxes would be useful for 
comparison with our slope estimates to determine potential 
evolutionary properties. We note again that the starburst 
and quiescent galaxy number-coun t slopes agree to within 
errors with the predictions made by Ranalli et al. (2003) us- 
in g both the X-ray/ radio correlation (e.g., Bauer et a l. 2002a' 
iRanalli et al.ll2003l) and the radio number counts of lRichardsl 

Visual extrapolations of the starburst, quiescent, and el- 
liptical galaxy number counts to fainter X-ray fluxes ap- 

We would prefer to use radio and infrared data, the "classical" probes of 
star-formation rate, to separate starburst and quiescent galaxies. However, the 
cun'ently available infrai'ed data cover only a tiny fraction of the CDF-N (e.g., 
[Alexander et al. 2002a), and regions with overlapp ing 1 .4 GHz and 8.5 GHz 
radio data (necessary for radio spectral slopes; e.g., Bauer et al. 2002a) either 
cover only a small fraction of the field (CDF-N) or are not yet deep enough 
to detect large numbers of star-forming galaxies (CDF-S). 



pear to be consistent with stacking results for various ex- 
tragalactic source populations (see Figure |9j. At the CDF 
flux limits, we find source densities of 984^[32, 343;|;g7, and 
105!^f sources deg-^ (soft) and 123!^;^, 48l!|^^, and 399![;5^ 
sources deg"^ (hard) for the starburst, quiescent, and elliptical 
galaxy samples, respectively. 

4.2. AGN Number Counts by X-ray Luminosity and 
Absorption 

Using our estimated redshifts, we can also divide the CDF 
sources into intrinsic X-ray luminosity and absorption bins. 
This is a useful exercise both for better understanding the 
types of X-ray sources contributing to the XRB a nd for test- 
ing the recent AGN X-ray luminosity function of lUeda et alJ 
(2003, hereafter simply the Ueda XLF). The Ueda XLF is 
constructed from X-ray sources with fluxes above F2-% kev = 
3.2 X 10"'^ erg cm"- s"' (including sources from the CDF-N) 
and is strictly applicable only for sources with intrinsic X- 
ray luminosities above L2-8 keV = lO**''^ erg s"' and column 
densities below A^h = lO^** cm"^. To compare to our full 
CDF dataset, we have extrapolated the Ueda XLF down to 
the CDF sensitivity limits and to intrinsic luminosities of 
1^2-% keV = 10"*°-^ erg s"' in order to test its predictive po wer. 
For details of this extrapolation, see lTreister et alJ J20041) . 

Figure^! shows the optimistic AGN sample in the CDFs, 
separated into several different intrinsic luminosity and ab- 
sorption classes. Also shown are the predicted X-ray number 
counts from the Ueda XLF for the same classes. Some cau- 
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Fi g. 10. — Plot identical to Fig.0except we show here a comparison of the soft (left) and hard (right) CDF number counts for all AGN (thick solid curves; 
see i|4. 1.21 and the predicted number counts from the Uedaetal. (2003) XLF (thin solid curves). The AGN number counts have been separated into several 
different intrinsic X-ray luminosity (top) and absorption (bottom) bins (color-coded following the order of the rainbow, with red denoting the highest luminosity 
and absorption bins) . 
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tion should be exercised when interpreting the CDF number 
counts fromA^H > 10^^ cm"^ AGN. Chandra's effective area 
above 4-5 keV declines rapidly such that the CDFs are almost 
certainly systematically missing a s ignificant fraction of the 
most highly obscured AGN (e.g., Treisteret al. 2004). This 
bias was not accounted for in our simulations, as we averaged 
over the energy dependence of Chandra's effective area in our 
flux correction and furthermore assumed that all of our sim- 
ulated sources had a spectral shape consistent with a F = 1 .4 
power law. 

It is readily apparent that the Ueda XLF overestimates the 
actual CDF AGN number counts in both bands. We have 
derived error bars on the cumulative distributions following 
lGehrelsl (|T986) and calculated the deviation at each data point 
from the model in units of sigma. In the soft band, the data 
points consistently lie below the model over the flux range 
2 X lO"'** erg cm"^ s"' to 5 x 10"'' erg cm"^ s"', with the sta- 
tistical deviation increasing from l-2cr at the brighter end of 
this range to 4-5a at the fainter end of this range. In the hard 
band, the data points consistently lie below the model over 
the entire flux range, with the statistical deviation steadily in- 
creasing with decreasing flux from 1-5.5(7. The discrepan- 
cies are still significant (up to 3cr) even if we only consider 
sources over the original flux range used to construct the Ueda 
XLF. Note that the upturn seen in the soft-band data below 
■fb.5-2.0 keV ~ 5 X 10"'^ erg cm"^ s"' in Figure lTUlis lik elv due 
to X-ray spectral complexity (as discussed in i]4.1.H not ac- 
counted for in the simple spectral model used to construct the 
Ueda XLF 

Turning to the individual intrinsic luminosity bins in Fig- 
ure ^1 we find that the soft-band XRB is dominated by 
sources with Lo.5-8 keV > 10'*^^ erg s"', while the hard-band 
XRB is dominated by somewhat less luminous sources with 
Lo.s-nkev = lO^^'^-lO^"''^ erg s"'. There is good agreement 
between the highest individual CDF luminosity bins and the 
Ueda XLF in both bands, with the best estimates in the 
soft band being among the Lo.5-8 keV > lO^'*''' erg s"' and 
Lo5-8keV = 10'*^'^-10'*^-^ erg s"' AGN samples, as weU as 
the bright end of the Lo.5-8 keV = lO'^^^-lO''^^ erg s"' AGN 
sample. The Ueda XLF, however, consistently overesti- 
mates the soft-band number counts of (1) Lo.5-8 keV = 10'*^'^- 
lO'*''-^ erg s"' AGN below 1 x 10"'"* erg cm"^ s~\ with the sta- 
tistical deviation steadily increasing with decreasing flux from 
1-6ct, and (2) Lo.5-8 keV < lO"*^-^ erg s"' AGN over the flux 
range 2 x 10"'^ erg cm"^ s"' to 5 x 10"'^ erg cm"^ s"', with 
the statistical deviation increasing from l-2fT at the bright- 
est and faintest ends of this range to 8-10(t in the middle of 
this range. In the hard band, the Ueda XLF again consistently 
overestimates the number counts of the (1) Lo.5-8 keV = 10'*^'^- 
2Q44.5 gj.g g-i ^Qjsj below 1 X lO"'"' erg cm"^ s"', with the sta- 
tistical deviation steadily increasing with decreasing flux from 
1-6(7, and (2) Lo.5-8 keV < lO'*^'^ erg s"' AGN over the flux 
range 2 x 10"'^ erg cm"^ s"' to 5 x 10"'^ erg cm"^ s"', with 
the statistical deviation increasing from 1-2(7 at the brightest 
and faintest ends of this range to 5-6(7 in the middle of this 
range. The discrepancies between the Ueda XLF predictions 
and the actual number of low-luminosity AGN have also been 
noted by Menci et al. (2004). 

In terms of intrinsic absorption bins, we find that in the soft 
band the XRB is dominated by relatively unobscured sources 
with A^H < 10^^ cm"^, while in the hard band the contribu- 
tion to the XRB from all of the different absorption bins is 
nearly equal. Again, there is good agreement between some 
of the individual CDF absorption bins and the Ueda XLF, 
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Fig. 1 1 . — We show a comparison of the distribution from spectral 
analysis of the CDFs (thin solid histogram) versus that predicted from the 
Ueda et al. 1 2003) XLF (thick solid histogram). Also shown are the 1 cr upper 
and lower limits of the CDF distribution (dashed and dotted histograms, 
respectively), indicating that the shape of the CDF distribution is relatively 
robust except for the lowest Wh bin. We caution that the actual distribution 
may be skewed towards higher Nw values if the redshift distribution we have 
assumed for the sources lacking redshifts is higher or our spectral models ai'e 
too simplistic (i.e., they fail to account for known spectral complexity in a 
systematic way). In the former case the shift should not be more than a factor 
of a few, while in the latter such spectral complexity would apply to both the 
CDF and Ueda distributions. 



although discrepancies between the Ueda XLF and the data 
are more evident here. In the soft band, the agreement is 
good for A^H < 10^' cm-2 and A^h = lO^^-lO^^ cm'^ AGN. 
However, the Ueda XLF model consistently overestimates the 
soft-band number counts of A^h = lO^'-lO^^ cm"^ AGN over 
the entire flux range, with the statistical deviation steadily 
increasing with decreasing flux from 1-19(7, and underesti- 
mates the soft-band number counts of A^h = lO^-'-lO^'* cm"^ 
AGN over the entire flux range, with the statistical deviation 
steadily increasing with decreasing flux from 1-6(7. Note that 
this latter discrepancy is in the opposite sense to Chandra's 
possible observational bias mentioned above and is there- 
fore likely to be even worse. In the hard band, the agree- 
ment between the Ueda XLF model and the data is good 
for A^H < 10^' cm-2, A^h = lO^^-lO^^ cm'^, and A^h = 10^^- 
10^"* cm"^ AGN, but consistently overestimates the hard-band 
number counts of A^h = lO^'-lO^^ cm"^ AGN over the flux 
range 2 x lO"''* erg cm"^ s"' to 4 x 10"'^ erg cm"^ s"' , with the 
statistical deviation steadily increasing with decreasing flux 
from 1-6(7. 

Some of the discrepancies between the individual luminos- 
ity and absorption bins and the Ueda XLF may be due to our 
assumptions regarding the redshift distribution (and hence the 
absorption and X-ray luminosity distributions). Several argu- 
ments suggest that this is unlikely to be the case: (1) The 
Ueda XLF overestimates the total number counts as well, 
where no redshift information is used. (2) The agreement be- 
tween the lowest and highest column density sources is fairly 
good. Significant changes to the redshift distribution either 
way would worsen the agreement. (3) If the redshift distribu- 
tion is on average higher than we have estimated in i|4.1l then 
sources from lower X-ray luminosity bins, where the over- 
estimate is already bad, would move into higher ones where 
the agreement is already acceptably good (thus worsening the 
agreement). (4) The X-ray sources are unlikely to lie at sig- 
nificantly lower redshifts than those assumed given that the 
sources are not generally detected in the bluest optical bands. 
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implying substantial redshifts or that they reside in severely 
underluminous galaxies. The latter is unlikely given what we 
know about the host galaxies of powerful AGN in the nearby 
Universe (Kauffmann et al. 2003). 

Figure shows the derived CDF A^h distribution (includ- 
ing errors) and the extrapolated Ueda et al. (.2003.) model. 
There is rough agreement between the data and model, al- 
though the CDF distribution clearly appears to be bimodel, 
while the Uedaetal. (2003) model does not. These devia- 
tions are less obvious in Figure [TOlbecause of the coarser A^h 
binning and the effect contamination from star-forming host 
galaxies has on the X-ray fluxes of highly obscured AGN. 
While the upper and lower limits on the overall A^h determina- 
tions (dashed and dotted lines, respectively) are tight enough 
to exclude the possibility of this bimodality being spurious, 
the assumed redshift distribution for sources lacking redshifts 
could play a role in shaping the A^h distribution. For instance, 
adopting a larger scatter in the redshift distribution might 
broaden the A^h distributi on to be mor e in line with the flat 
distribution from the Ueda et al. (2003) model. Alternatively, 
if we have systematically underestimated redshifts this would 
tend to shift the second A^h peak toward higher values and 
increase discrepancies between the data and model at mod- 
erate column densities. We consider the CDF A^h distribution 
shown in FigurelTTIto be fairly robust, however, since it would 
take exceptionally large deviations in the adopted redshift dis- 
tribution to change intrinsic A^h values by more than a factor of 
a few given that Nwiintr) « ( 1 +1)^ -INniobs). Thus, it appears 
that the extrapolated Ueda et al. (2003) A^h distribution is not 
entirely appropriate for the CDF data. Treister et al. ("2004*) 
come to a similar conclusion using a more luminous subset 
of the CDF AGN and an A^h distribution derived instead from 
X-ray colors. 

4.3. Contributions to the Extragalactic X-ray Background 

The contributions to the soft and hard XRB from each 
CDF source type are shown in Table |2l normalized by the 
average total flux density of the XRB [taken from M03 in 
the soft ba nd to be (7.52 ±0.35) x 1 0"'^ erg cm"^ s"' deg"^ 
and from De Luca & MolendU 120041 in the hard band to be 
(2.24 ± 0.11) X 10"" erg cm"^ s"' deg"^]. For sources 
brighter than observed in the CDFs (i.e., Fo.5-2.0 kev > 8 x 
10"'"^ erg cm"^ s"' and fa-g keV > 1 x 10"'-^ erg cm"^ s"'), we 
have adopted total resolved flux-density contributions from 
the models of M03, which are 19.2% in the soft band and 
1 1.5% in the hard band; we assume uncertainties of ^4% on 
these values estimated from Figure 5 of M03. 

The crude breakdown of these bright resolved fractions is 
40% stars, 40% Type 1 AGN, 5% Type 2 AGN, 10% clus- 
ters, 5% galaxies in the soft band, and 1% stars, 69% Type 1 
AGN, 20% Type 2 AGN, 8% clusters, and 2% galaxies in the 
hard band (e.g.. iKrautter et al . 1999; L a Franca et alJl2002t 
lAkivama et alJl2003l) . When known, these values have been 
added to column 4 of Table |2]to produce total resolved XRB 
fractions in column 5 of Table |2l In total, and neglecting 
possible large-scale structure effects and instrumental cross- 
calibration uncertainties, we find 89.5:l:g 7% and 86.9:l:g f % of 
the XRB has been resolved in the soft and hard bands, respec- 
tively. 

We find that AGN as a whole contribute «83% and «95% 
to the resolved soft and hard XRBs, respectively. Dividing the 
AGN further, we find that Type 1 AGN alone contribute sig- 
nificantly (~ 57%) to the AGN fraction of the resolved XRB 
in the soft band, but far less so ('^ 35%) in the hard band. X- 



ray unobscured/mildly obscured AGN completely dominate 
('- 76%) the AGN fraction of the resolved XRB in the soft 
band, and contribute nearly half (~ 46%) in the hard band. 
If this latter trend continues to harder energies, it suggests 
that highly obscured X-ray sources will likely dominate the 
XRB emission at higher energies (e.g., closer to w 20^0 keV 
peak of the XRB). This can be directly tested by perform- 
ing similar analyses on deep 5.0-10.0 keV number counts 
from XMM-Newton. In contrast, star-forming galaxies com- 
prise only ?alO% and «2% of the soft and hard XRBs. Star- 
bursts appear to make up the bulk of the X-ray emission in 
the soft band, while quiescent galaxies apparently dominate 
in the hard band; contamination from embedded, highly ob- 
scured AGN may be responsible in the latter case. 

Extrapolating the total number-count slopes from Table ^ 
down to 1 X 10"''^ erg cm"^ s"' in the soft band and 1 x 
10"'^ erg cm"^ s"' in the hard band can account for an ad- 
ditional 6.6!^ j% and 24.l!^|^% of the XRBs, respectively. 
We note, however, that some of the number-count slopes for 
the individual CDF source types are significantly steeper than 
the overall number-count slopes and that independent extrap- 
olation of these slopes actually exceeds the total XRB by a 
large percentage in some cases. This suggests that (1) known 
point-source populations can explain all of the remaining un- 
resolved XRB and (2) the steep slopes found for some source 
types must flatten in the decade of flux below the CDF sen- 
sitivity limits. As such, we contend that there is little room 
for extra contributions from, e.g., a truly diffuse compo- 
nent. We additionally find that the number density of star- 
forming galaxies should overtake that of the AGN just below 
F{).5-2.Q keV ~ 1 X 10"'^ erg cm"^ s"' in both bands, and these 
sources should contribute an additional ^10% and ^5% to the 
soft and hard XRB emission (here we have assumed a flatten- 
ing of the starburst number counts fo r sources as inferre d from 
infrared galaxy number counts; e.g. JCharv & Elbaz"2001). 

Thus, the overall contribution of galaxies to the XRB ap- 
pears to be 10-20% in the soft band and 2-6% in the hard 
band, indicating that if the cosmic star-formation rate evolves 
as ( 1 + zy then q is likely to be near q = 3 rather than a 
higher value (see Persic & Rephaeli '2003"). This is con- 
sistent with previous results from stacking analyses (e.g., 
iHornschemeier e t al. 2002; Ge orgakakis et al . 2004) a s well 
as the recently measured normal galaxy XLF CNorman et alJ 

We note that we are using relatively broad X-ray bands, and 
that the resolved fraction at a given energy may change some- 
what throughout these bandpasses. Since the effective area 
of Chandra falls off rapidly above 4 keV, we are likely to 
resolve more of the 2-4 keV XRB and less of the 4-8 keV 
XRB. An investigation into this is beyond the scope of this 
study, but has been undertaken by Worslev et al. (2004; 200^ 
in preparation) using the deep Lockman Hole and the CDFs. 
They perform photometry of the resolved source in several 
narrow X-ray bands, finding that the resolved fraction of the 
XRB is indeed sytematically lower at harder X-ray energies. 
They also find that the resultant X-ray spectrum of the un- 
resolved XRB is consistent with that from a highly-obscured 
(A^H ~ lO^^-lO^"* cm-2 for z < 2) AGN. 

5. CONCLUSIONS 

We have combined the CDF-N and CDF-S X-ray sam- 
ples, along with substantial public ancillary data, to exam- 
ine the X-ray number counts as a function of source type. 
Extensive simulations were carried out to quantify and cor- 
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rect for the completeness and flux bias problems which affect 
the CDF number counts. Once corrected, we find that the 
number counts from the two fields are consistent with each 
other aside from sources detected in the 2-8 keV band below 
■P2-8 keV ~ 1 X 10"''' erg cm"^ s"', where statistical deviations 
gradually increase to 3.9a at the faintest flux levels. We also 
find that our overall number counts are consistent with previ- 
ous determinations. In total, we have resolved ^9.5'!ifj% and 
86.9!^:^% of the extragalactic 0.5-2.0 keV and 2-8 keV XRB, 
respectively. 

Using a classification scheme based on X-ray spectral 
properties, intrinsic X-ray luminosities, radio morphologies, 
variability, optical spectroscopic classifications, and X-ray- 
to-optical flux ratios, we have separated the CDF X-ray 
sources into 698 AGN (split into optical Type 1 and not ob- 
viously Type 1, and X-ray obscured and unobscured/mildly 
obscured), 109 star-forming galaxies (split into starburst, qui- 
escent, and elliptical), and 22 Galactic stars, and determined 
their individual number counts. We additionally calculate 
the number-count slopes and normalizations below 2 x 10"'^ 
erg cm~^ s~' for all source types assuming a single power-law 
model. 

We confirm that AGN power the bulk of the extragalactic 
XRB, with only a small contribution from star-forming galax- 
ies. The most significant contributions to the XRB are from 
sources with Lo.5-8 keV > lO'*-''^ erg s"' and A^h < 10^~ cm"^ 
in the soft band and Lo.5-8 keV = lO'^^^-lO^"^^ erg s"' and an 
evenly-distributed range of absorption column densities in 
the hard band. This trend suggests that even less luminous, 
more highly obscured AGN may in fact dominate the num- 
ber counts at higher energies, where the XRB intensity peaks. 
At the CDF flux limits, the overall AGN source densities are 
7166;|;2*92 sources deg~^ and 4558;I;2()7 sources deg"^, respec- 
tively, factors of ^ 10-20 higher than found in the deepest 
optical spectroscopic surveys. 

While star-forming galaxies make up a small fraction of 
sources with fluxes higher than ^ 10"'^ erg cm"^ s"' in 
both bands, their numbers climb steeply below this flux 
such that they eventually achieve source densities of 1727;|;j^g 
sources deg"^ and 71 lilij,™ sources deg"^ at the CDF flux lim- 
its (with starburst galaxies making the largest contribution) 
and comprise up to ^ 40% of the sources at the faintest X-ray 
fluxes. Extrapolation of the number-count slopes for galaxy 
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source types can account for all of the remaining unresolved 
soft and hard XRBs. Moreover, within a factor of a few be- 
low the current CDF soft-band flux limit, the sky density of 
star-forming galaxies will likel y overtake that of A GN. 

We have also compared the Ueda et al. (2003) XLF with 
our X-ray number counts. While the agreement is generally 
good for most subsets of CDF sources, it appears that the ex- 
trapolated XLF requires some significant refinements before 
it will be able to reproduce the data with sufficient accuracy. 

Further improvements to our study can be made in several 
ways. Constraining the bright-flux end of the number counts 
for all of our source types would be particularly useful for de- 
termining more accurate slopes and XRB contributions. Ad- 
ditional observations of the 2 Ms CDF-N with Chandra would 
(1) yield additional photons to discriminate better between 
AGN and galaxies and model spectral complexities for bet- 
ter absorption column density and luminosity estimates, (2) 
allow confirmation of the sharp rise in the number counts of 
both star-forming galaxies and moderate-to-high obscuration 
(A^H 10^2- lO^"*) AGN that we find here, and (3) perhaps 
provide evidence for their eventual flattening as is required 
by the overall XRB flux density. Upcoming ancillary obser- 
vations at infrared (Spitzer) and radio (VLA) wavelengths, as 
well as detailed analyses of the existing optical spectra, would 
be extremely helpful for refining our AGN and galaxy sub- 
classifications. Finally, a high-sensitivity, high spatial resolu- 
tion X-ray telescope able to probe to 10"'^' erg cm"^ s"' in the 
10^0 keV band is desperately needed to resolve the peak of 
the XRB and determine once and for all the nature and com- 
position of the XRB. 
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